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Abstract 

For neutrino oscillations to take place the entangled quantum state of a neutrino and a charged 
lepton produced via charged current interactions must be disentangled. Implementing a non- 
perturbative Wigner-Weisskopf method we obtain the correct entangled quantum state of neutri- 
nos and charged leptons from the (two-body) decay of a parent particle. The source lifetime and 
disentanglement length scale lead to a suppression of the oscillation probabilities in short-baseline 
experiments. The suppression is determined by TxLsjLoac where is the smallest of the decay 
length of the parent particle or the disentanglement length scale. For Lg > Lose coherence and 
oscillations are suppressed. These effects are more prominent in short base line experiments and at 
low neutrino energy. We obtain the corrections to the appearance and disappearance probabilities 
modified by both the lifetime of the source and the disentanglement scale and discuss their impli- 
cations for accelerator and reactor experiments. These effects imply that fits to the experimental 
data based on the usual quantum mechanical formulation underestimate sir? [29) and (5m^, and are 
more dramatic for <5m^ ~ eV^, the mass range for new generations of sterile neutrinos that could 
explain the short-baseline anomalies and long disentanglement length scales. 

PACS numbers: 14.60.Pq;13.15.+g;14.60.St 
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I. INTRODUCTION 



Neutrino masses, mixing and oscillations are the clearest evidence yet of physics beyond 
the standard model [l|-[3[. They provide an explanation of the solar neutrino problem [^|6[ 
and have important phenomenological (iMsl. [?- astrophysical j5,|l2, 13 1 and cosmological 



14l | consequences. A remarkable series of experiments have confirmed mixing and oscillations 
among three "active" neutrinos with 5m^ = 10~^ — 10~^eV^ for atmospheric and solar 



oscillations respectively [151]. 

A fascinating aspect of neutrino oscillations is that they provide an extraordinary example 
of macroscopic quantum coherence maintained over hundreds of kilometers. It is particularly 
this aspect that has sparked an ongoing discussion in the field that seeks to clarify the main 
concepts behind the physical interpretation of oscillations. 

As neutrino oscillations open a window to explore physics beyond the standard model, it 
is important to understand the underlying phenomena at the deepest level, and the domain 
of validity of the various calculations of oscillation probabilities and their impact on experi- 
ments. In particular, the standard approach of treating neutrino oscillations by analogy with 
Rabi-oscillations in a two state system (see for example, [l|-[l, 0, [sl and references therein) 
while simple and intuitive, has motivated a wide ranging discussion. Deeper investigations 
of this basic paradigm have already raised a number of important and fundamental questions 
lM3 that are still being debated |20|. 



A correct interpretation of the results from oscillation experiments require understanding 
of both the production and detection mechanisms (21 - 23| . Neutrino detection is indirect 
through charged or neutral current processes and mostly through the detection of the as- 
sociated charged lepton. As for the production mechanism, the neutrino state is produced 
by the decay of a parent particle via charged current interactions. Coherence (and deco- 



herence) aspects of the production and detection of neutrinos |23l . l25l . |26| and lifetime of 
the source |2J] have been discussed, however only recently the recognition that the neutrino 
state produced by the decay of a parent particle via charged current interactions is in fact 
entangled with that of the charged leptonjisl. 27-33] has become the focus of a reassessment 
of the dynamics of neutrino oscillations. 

Quantum entanglement is a direct consequence of conservation laws in the production 



process [18|, l27j which result in a correlated quantum state of the neutrino and its charged 
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current lepton partner. As observed in refs. [2^ 
to interfere coherently and oscillate, the quantum state must be disentangled: entanglement 
surviving for a very long time projects out states of definite energy and prevent oscillations. 
In a typical experiment disentanglement of the charged lepton occurs when this particle is 
measured, absorbed or decays, after disentanglement the quantum state is reduced and is re- 
set. The full dynamics of the process of production of the entangled state, disentanglement 
and further evolution to production of another charged lepton at a (far) detector was studied 
in a model in ref.jl^ with a focus on long baseline experiments. Ref.jl^ studied the free 
time evolution of a disentangled wave-packet produced from pion decay including lifetime 
effects but without addressing the production of charged leptons by the disentangled state 
and detection process . 

The results of ref. show that if disentanglement occurs on time scales much shorter 
than the oscillation scale, namely for long baseline experiments, the familiar result obtained 
from the simple quantum mechanical picture and factorization is reproduced, but also point 
out possible subtle consequences if the disentanglement process occurs on time scales of 
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the same order of or longer than the oscillation time scale, with potential impact on short 
baseline experiments. This possibility is also hinted at in ref. [31|. 

In the last few years several experimental results have been accumulating that cannot be 
interpreted within the "standard paradigm" of mixing and oscillations among three "active" 



neutrinos with ~ 10~^ — 10^'^. The early results from the LSND experiment 351 ] have 
recently been confirmed by MiniBooNE running in antineutrino mode[36| both suggesting 
the possibility of new "sterile" neutrinos with 6m'^ ~ eV^. More recently, a re-examination 
of the antineutrino flux ^3] in anticipation of the Double Chooz reactor experiment resulted 
in a small increase in the flux of about 3.5% for reactor experiments leading to a larger deficit 
of 5.7% suggesting a reactor anomaly\^^. If this deficit is the result of neutrino mixing and 
oscillation with baselines L < 10 — 100 m, it requires the existence of at least one sterile neu- 



trino with (5m^ > 1.5 eV^ and mixing amplitude sin^(26') ~ 0.115[38[. Taken together these 
results may be explained by models that incorporate one or more sterile neutrinos that mix 



with the active ones |39H42| including perhaps non-standard interactions |43|. Furthermore 
the latest analysis of the cosmic microwave background anisotropics by WMAPji^] suggests 
that the effective number of neutrino species is N^ff = 4.34 ± 0.86 and XK^f) < 0.58 el^ 
bolstering the case for sterile neutrino(s) with m < eV. 

The common aspect of accelerator and reactor anomalies is that these are all short baseline 
experiments with 10m < L < 600m and this aspect, when considered along with the 
potentially relevant corrections from disentanglement discussed above, motivate our study 
of the disentanglement and lifetime effects on short baseline oscillations. Although the effect 
of the muon lifetime on neutrino oscillations has been studied in ref. j45| and more recently 



ref.[46[ argued that the pion lifetime introduces decoherence on oscillations — at the 
near detector of the MINOS experiment, the combined effect of lifetime of the source and 
disentanglement of charged leptons have not yet been discussed with regard to the distortion 
of the spectrum for charged lepton events at the detector in short baseline experiments. 

Goals: The accumulation of experimental evidence of short baseline anomalies and the 
recognition that entanglement and lifetime effects may lead to corrections in the oscillation 
probabilities precisely in short baseline experiments motivates us to understand the im- 
pact of these effects on the appearance and disappearance probabilities, and their possible 
experimental implications. 

In this article we seek to understand the subtle aspects arising from quantum mechanical 
correlations as a result of the fact that the neutrino states produced at charged current 
vertices are entangled with the charged lepton partner. Measurement, absorption or decay 
of this charged lepton disentangles the quantum state, but the emerging neutrino state 
carries information on the quantum correlations in its evolution. These correlations along 
with intrinsic energy uncertainties associated with the lifetime of the parent particle whose 
decay produces the neutrinos, influence the oscillation probabilities. Our goal is to study 
these corrections in the simplest and most clear setting that allow a systematic calculation of 
the effects and to extract the possible impact of these effects on the experimental observables 
and their interpretation. 

Results: 

The dimensionless parameter that determines the impact of the lifetime of the source and 
the disentanglement time scale on the oscillation probabilities is 

TxLs/ Lose 
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where Ls is the smaller between the decay length of the source and the disentanglement 
length scale at which the charged lepton produced with the neutrino is measured, absorbed 
or decays and Lose oc Ej,/5rn? is the oscillation length. 

A detailed analysis of the production and disentanglement of the quantum states of 
neutrinos produced in charged current vertices reveals that the usual "Pontccorvo" states 
familiar from the simple quantum mechanical approach, are a reliable description only when 
Lg <C Lose in which case appearance and disappearance probabilities are given by the usual 
expressions, but for Ls 3> Lose the energy uncertainties become smaller than the distance 
between energy eigenstates and coherence is suppressed with a concomitant suppression of 
the appearance probability. 

Appearance probabilities are suppressed in short baseline experiments where both the 
lifetime of the source and the disentanglement scale are comparable to or a large fraction of 
the baseline. 

We find that at MiniBooNE the dominant source of suppression is the decay length of 
the pions, whereas at LSND and reactor experiments we argue that the relevant length scale 
is the disentanglement distance. 

In all these cases, fits of the experimental data to the usual quantum mechanical appear- 
ance and disappearance probabilities underestimate both dm? and sin^(2^). These fits are 
much less reliable at low neutrino energy (for fixed Srn^ , L) because for low energy events 
the ratio Lg/ Lose is larger and the suppression of the oscillation probability is stronger. 

The corrections from lifetime and disentanglement efi'ects on short baseline experiments 
are more dramatic for the mass range (5m^ ~ 1 eV^ which is the putative mass range for 
sterile neutrinos that could explain the short baseline anomalies, and for low neutrino energy. 



II. A MODEL OF "NEUTRINO" OSCILLATIONS 

In order to exhibit the main results in a clear and simple manner, we introduce a bosonic 
model that describes mixing, oscillations and charged current weak interactions reliably. 
The complications associated with fermionic and gauge fields are irrelevant to the physics 
of mixing and oscillations, as is obviously manifest in the case of meson mixing. 

We study the model defined by the Lagrangian density 

£ = >Co[W^, ^a] +>CoM +Ant[7r, Za,i^a] ; « = 6, /X (II.l) 

with 

Cq[u] = - [a^^^a^^ - ^^M^] , (II.2) 

where ^ is a fiavor doublet representing the neutrinos 

and M is the mass matrix 



(II.4) 



The interaction Lagrangian is similar to the charged current interaction of the standard 
model but explicitly includes a vertex that describes the decay of a parent particle (here the 
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pion) into a charged lepton and its flavor neutrino vr — )■ /i z/^, namely 



i^\nt{x, t) = 5'^7r(f, t)/i(f , t)u^{x, t) + g W{x, t) e{x, t) z/e(x, t) + ^{x, t) z/^(x, t) , (II.5) 

where g plays the role of the electroweak charged current coupling, and gT, includes the 
pion decay constant. 7r(x) represents the pion fleld or alternatively any parent particle that 
decays into a charged lepton and its associated neutrino, W{x) represents the vector boson, 
and I = e, the two charged leptons. 

Obviously this simple model cannot describe CP violating effects or distinguish between 
neutrino vs. antineutrino modes, however our goal is to understand decoherence effects 
associated with lifetime of the source and disentanglement of charged leptons. 

The mass matrix is diagonalized by a unitary transformation 

u-Wio) . [-^ I) ; uio) . { ^-^^ Z'e) ■ ("■«) 

In terms of the doublet of mass eigenstates, the flavor doublet can be expressed as 

u{e) (11.7) 



This bosonic model clearly describes charged current weak interactions reliably as it includes 
all the relevant aspects of mixing and oscillations. Furthermore the coupling vr /i z/^ allows 
us to study the dynamics of the production process including the lifetime of the source ( 
pion) within the same model. 

In order to clearly separate the effects from the lifetime of the source and disentanglement 
time scale from the effects of wave-packet localization, this article is primarily devoted to 
the analysis in terms of plane waves, in section (jV]) we comment on the modiflcations from a 
wave packet treatment, but postpone the full treatment with wavepackets to a more thorough 
forthcoming study [13]. 

We consider the case in which a neutrino and its flavor charged lepton partner are pro- 
duced via the decay of a parent particle, in this case a pion, however, the discussion and the 
main consequences are general, with the only difference being the associated many-particle 
phase space if the decay is in more than two particles. The production process corresponds 
to TT — ?► yU z/^ where the fi is "observed" or is absorbed (or decays) at a "disentanglement" time 
scale tfj^, the disentangled neutrino is detected via a charged current interaction ^ Wl 
where / = e, /i is the charged lepton. 

The Wiener- Weisskopf |5 71. method described in the appendix yields the entangled 



state that results from pion decay, the relevant part of the interaction Hamiltonian in the 
interaction picture is 

H^[t) = 9^ j 7r(f , t)ii{x, t)Uf,{x, t) (II.8) 

where the time evolution is that of free flelds. The neutrino field operator 

z/^(x, t) = cos{9)h'2{x,t) — sm{9)h'i{x,t) , (II-9) 

where i'i^2{x, t) are expanded, as usual, in annihilation and creation operators of mass eigen- 
states ai 2(p); a\ 2{p) respectively, with the single particle mass eigenstates being 

4im = Wp.) (11.10) 
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where |0) is the vacuum state annihilated by ai{p). We note that the transformation law 
(]II.7|) applies to the field operators not to the single particle states. 

Consider that the initial state, at t = is given by a single particle pion state described 
by a plane wave with momentum k, namely 



1^(^ = 0)) = kg). 

Hf connects the single particle pion state to the states |k) 
matrix element is given by 

^-i{E^{k)-Ef,iq)-Qiip))t 



(/Xg-, i^j;,i\HJ{t)\U^) = Mi{k, q,p, t) 



V 



(11.11) 

fJ'q,)Wp,i)- The transition 



(11.12) 



where fli{p) = a/p^ + mf are the energies of the neutrino mass eigenstates. 

We are now in position to use the results of the appendix for the time evolved state 
in the Schroedinger picture resulting from pion decay, it is given by eqn. OA. 271) with 
Cyi(O) = C^(0) = 1 and the set k described above, we find 



where the amplitudes 

Ci{k,q,p,t) = Mi{k,q,p,t = 0) 



5^C.(fc,g,p,t)e-^(^''(^)+^'(^))*|/i,--) 1^.,^,) ; q = k-p, 

(11.13) 



Es~ni{p)-^T^{k) 



Es = E^{k)-E^{q) (11.14) 



E^^lk) is the fully renormalized pion energy (including the self-energy correction from the 
intermediate states, see appendix) and r^(fc) = M^ro/ii^^(fc) where To is the decay rate of 
the pion at rest. For = the entangled state fill. 131) is the same as that obtained in 
ref. 291] in lowest order in perturbation theory. 

Although the particular form of Aii{k,q,p,t) for the bosonic theory considered here is 
given by (111.121) the results ( 1II.13|II.14I) are general in terms of the transition matrix element 
Miik,q,p,t = 0). 

The state (]II.13|) is an entangled state of the neutrino mass eigenstates and the muon, the 
entanglement is evident in that it is a sum of product states, not a simple product state, the 
amplitudes C are a measure of the correlation between i/j and the charged lepton. The en- 
tanglement is a consequence of momentum conservation since q = k—p. The "observation", 
measurement or decay of the muon state at time t^ disentangles the neutrino state. If the 
muon is "measured" in a plane wave state with momentum Q the disentangled state is ob- 
tained by projecting the quantum state fill. 131) onto the state l/i^), namely, the disentangled 
neutrino state is given by 



|v^(t^)) = 



-iE^(Q)t^ 



WE^{k)E^{Q) 



j=l,2 



Fi[k,Q,p;t^]\v^^i) ; p = k-Q 



where 



Fi[k,Q,p;t^ 



I _ g-i{Ss-ni{p))tp g-r^{fc)tp/2 
Es - a(p) - iK{k) 



^5 = E^ik) - E,,{Q) 



(11.15) 
(11.16) 
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The functions Fi encode the information of production of the entangled charged-lepton- 
neutrino pair and the measurement of the charged lepton, it features both time scales: the 
lifetime of the source and the disentanglement time scale t^. 

The number of muons of momentum Q detected at is given by 



dVQ 



{^{t,)\al{Q)a,{Q)\^{t,)) = {VM\VM) , (11.17) 



where a|j, are creation and annihilation operators for muons. Thus the normalization of 
the disentangled neutrino state is completely determined by the number density of muons 
detected at t^. 
We find 



9l 



4VE,{k)E^{Q) 



F2{k,Q,p,t,)\' + 



F,{k,Q,p,t^)f 



This expression becomes familiar from the following analysis: the functions 



\Fi{k,Q,p,t^)[ 



1 — e 2*Mj+4e sm 



(^5 - n,{p)) 



r2(fc) 



(11.18) 



(11.19) 



are strongly peaked at Es = Ej^i^k) — -E^(Q) = ^i{p) with width determined by the largest 
of TT^{k) ; 27r/t^ becoming proportional to energy conserving delta functions in the limit 
when these become very small. This is clearly seen in two relevant limits: 



i): the narrow width limit with Tj^lk) <C 1 but large where 

„ Asm'[{Es-niip)yfj 



\Fi{k,Q,p,t^)[ 



(^Es-n,{p))y 



(11.20) 



This limit corresponds to a long disentanglement time scale but <^ T^^ik) where TT^{k) = 
l/r^(A;) is the pion lifetime in the laboratory frame. This function is strongly peaked at 
Es = ^i{p) with maximum height and width 2T{/t^ which is the largest of VT^{k) and 
27r/t^ for this case. As t^ becomes very large, 



\Fi{k,Q,p,t^)\^ ~ 2Tit^5{Es - 0,(p)) . 



(11.21) 



ii:) the opposite limit, for VT^{k)t^ ^ 1, in which the disentanglement time scale is 
much longer than T-^{k), where 



|F,(A;,Q,p,i^)f~ 



Es-^iip))' 



+ 



(11.22) 



The function is strongly peaked at £"5 = VLi{p) of height A/V\{k) and width r^(A;)/2 which 
is the largest of ; 27r/i^ in this case. In the narrow width limit 

|2 



\Fi{k,Q,p, tf. 



27rT^{k)5{Es-ni{p)). 
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(11.23) 



It is clear from this discussion that describe approximate energy conservation at the 
production vertex, approximate because the finite disentanglement time scale and/or the 
pion lifetime T.,r{k) broaden the energy conserving delta functions with an energy resolution 
determined by the width which is the largest of 2n/t^ or 27r/Tn-(/i;) respectively, namely the 
shortest time scale. 

For the general form (III.lQp . a straightforward integration yields 



iFA^dEs 



27TT^{k) 



(11.24) 



therefore, assuming that the energy distribution is very sharply peaked at Eg = fij(p), we 
can approximate 



\Fi{k,Q,p, t^ 



2nT^{k) 



1-e 



6{Es-n,{p)) 



(11.25) 



In this approximation the total number of muons measured at the disentanglement time 
scale is 



N„ 



(fx 



(27r)3 



Af^iQ,t^)c^T^{k)T^ik) 



where 



r^(A;) = cos^(6') T^^^^^{k) + sin^(6')r^^^^j(A;) 
is the total pion decay rate (in this simple model) and 

d^Q 



(11.26) 
(11.27) 



(k) 
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32 7:^ E^{k) J E^{Q)ni{p) 



5{E^{k)-E^{Q)-n,{p)) ■ p = k~Q (11.28) 



are the partial widths for pion decay into the neutrino mass eigenstates. Although this result 
applies to the simple model considered here, clearly it is general and conceptually correct: 
the particle decays into mass eigenstates which are the correct eigenstates of the unperturbed 
Hamiltonian, with the probabilities determined by cos^(6') , sin^(^) respectively. 

In what follows, we consider ultrarelativistic and nearly degenerate neutrinos and write 



Qiip) - E,{p) - A{p) ■ n2{p)^E,{p) + A{p), 



where 



Eu{p) 



o m? + rrin 
P + ^-1^ — - 



A{p) 



Srn^ 



AE,{p) 



6m 



m 



1) 



(11.29) 



(11.30) 



and take A{p) ^ Ei,{p) as is the experimentally relevant case. For AjEy ^ 1 the energy 
conserving 6{E^{k)-E^{Q)-^i{p)) in fllL28|) maybe replaced by 5(E^(A;) -E^(Q) -E^(p)) 
in the integral, because for the experimental range bm^ < 1 eV^ ; Ey>\ MeV the relative 
error incurred oc [A/Ey) < 10~^^ is much smaller than typical experimental resolution. 

The detection or measurement of the muon at time re-sets the quantum state to 
|V^(tu)), upon further evolution in time this disentangled state evolves into 



-iHot 



U{t,t^) 



(11.31) 
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where 



(11.32) 



is the time evolution operator in the interaction picture with boundary condition W(t^, t^) = 
1. 

The usual "Pontecorvo" quantum state familiar in the literature are simple linear super- 
positions lui) = Yli^iiWi)^ where Uu are the elements of the mixing matrix OII.Gp . and the 
corresponding muon neutrino state at any time t is 



|z/^)(t) = cos(^) e-^^^^*')* |z/2) - sin(^) e'^^^^^)* . 



(11.33) 



Instead, the corresponding quantum state evolved freely in time from up to t from the 
disentangled state |V^(tu)), is obtained from fill. 311) by setting W(t,t^) = 1, it is given by 



|V^(t)) = e-^^"(*-*'^)|V^(t^)) 



N 



cos{9) F2{k,Q,p,t^) e-^^^(P)* 



- sm{9) F^{k,Q,p,t,) e-^^^^^')* |z/p,i) 



(11.34) 



where the pref actor can be read off (III.ISP and we have neglected terms of 0(A(p) / E^{p)Y 
thereby approximating a/2 ^i{p) — a/2 E^ij)) in fill. 151) including this factor in A^. Conse- 
quently we also replace VLi^2{p) ~^ Ey{p) in flll.isp . 

Obviously the usual "Pontecorvo" states emerge up to the overall normalization factor 
if \Fi\ = 1 since time independent phases can be absorbed in the definition of the mass 
eigenstates. The conditions under which this equality is fulfilled is analyzed below. 

Following the familiar quantum mechanical approach to obtain the survival probability 
we find 



'Pu,^u,{t;t,) = \{V,{t,)\V,{t))\ 



cos 



F2P + sin^ 



^il 



- |F2nFi|2 sin" (29) sin^ 



4^.(p) 



(t-t,) (11.35) 



Obviously if = IF2I there is agreement with the usual result from Pontecorvo states 
up to an overall normalization. 

We note that invoking the approximation (]II.25|) for yields IF1F2P = as the 
product of delta functions vanishes for Qi 7^ Q2 thereby leading to the hasty conclusion that 
coherence is completely suppressed, however (]II.25P is an approximation that neglects the 
fact that are not sharp distributions but broadened with typical widths F^r or 27r/t^. 
Thus an assessment of the coherence leading to oscillations and interference requires a careful 
and detailed examination of the product |FiF2p. 

It proves convenient to use (1II.16|II.29P and write 



\Fi{k,Q,p;t^) 



r7r(fc) , 

e 



rv(fc) . . n 

+ 4 e 2 sm 



(Ss + Aip))'-^ 



(^5 + A(p))' + 



r|(fe) 



(11.36) 
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fl_e-^'A +4 6-^*" sin^ \{£s-A{p)yf 

\F2{k,Q,p;t^)\' = ^ ^- ^ (11.37) 

riik) 



(Ss-Aip)^ + 



where 

Ss = Es-E,{p). (11.38) 
The hmits studied above clarify the impact of the width of the energy distribution, 



i): T.„{k)tfj_ ^ 1. In this case 



r,,(fc) 

,2 V 2 



F,{k,Q,p;t^)\ ^2T^{k) ^ -{ (11.39) 



r|(fc) 



\F,{k,Q,p;tX -'^T^ik) ^ (11.40) 



£s + A{p) 



£s-Aip)] +^ 



where T^^lk) = l/r^(A;) is the pion hfetime in the laboratory frame. Each Lorentzian 
is peaked at Ss = ±^(p) respectively with a height oc T^{k) and width oc r^(/c), their 
product is depicted in fig. ([1]) for T.„{k) > A(p) and r^(fc) ^ A(p) respectively. For 
r^j-lk) ^ A(p) the product is similar to one Lorentzian peaked at £^5 ~ because the 
width of the individual Lorentzians (~ rT^{k)) is larger than their separation (~ A(p)), 
therefore the lifetime of the source introduces a large energy uncertainty that cannot 
resolve between the nearly degenerate energy eigenstates and blurs the individual 
peaks under one broad peak. In this case ~ and the disentangled state 

|V^(t)) is proportional to the corresponding Pontecorvo state. 

On the other hand if r^(A;) ^ A{p) the product is a double peaked distribution with 

peaks at Ss ^ ±A(p) of widths ~ r^(A;) but with height 1/ (^A'^{p)rl{k)^ < '^/KikY 

as is the case for |-Fi,2|^ which are the terms that do not feature oscillations. In this 
case the distance between the peaks ~ A{p) is much larger than the width of the 
individual peaks ~ r^(/c) and the energy uncertainty ~ r^(/c) is small enough that 
the individual mass eigenstates are resolved. 

This can be seen more efficiently from the identity 



£| + A» + 3^ 2r^'' + '^^n "2A2(p) + ^ 2 



1 1 



2 



|2 



(11.41) 

where the last approximate equality follows from the fact that |-Fi,2p are strongly 
peaked at Ss T A{p). 

Therefore, the conclusion is that when A{p) 3> r^(/c) and 3> it follows that 
|Fip|F2p -C 1-^1,21^ and coherence and oscillations are suppressed. 



^ This suppression is not manifest in fig. ([T]) because we have chosen A as the overall scale for presentation 
purposes. 
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ii): r^(A;)t^ ^ 1. In this case 



\Fi{k,Q,p]t^)\ 



\F2{k,Q,p;t^ 



2 ''A' 



+ 4 sin^ 



(£5 + A(p))| 



4 



r^(fc) , 

2 ''A' 



+ 4 sin^ 



{Ss-A{p)) + 



r|(fc) 



(11.42) 



(11.43) 



These are sharply peaked distributions that feature maxima at £^ = ^A{p) respec- 
tively with heights and widths ~ '^'n'/tfj,, displayed in fig. ([2]) for A{p) = 1; T.„{k) = 
0.01; tfj_ = 2,20. If 27i/t^ > A{p) the two peaks are blurred into one broad peak, 
whereas if 27r/t^ ^ A{p) the peaks are separated, this is a manifestation of the sepa- 
ration of mass eigenstates in real time as a consequence of the energy-time uncertainty, 
nevertheless, the product |FiF2p is not only not vanishing but with large support at 
the individual peaks. 

An important case is when r,r = 0, namely a stationary source, which will be important 
in the discussion below. In this case for large 



sm 



1^1, 



{Ss±A{p)y- 



Ss±A{p) 



27it^6{£s±A{p)) 



(11.44) 



whereas 



1-^1-^2 



sin[(g5-A)|] smjiSs + AYf] 
i£s - A) 
sin{A tf. 



i£s + A) 



+ 2 



A 

cos(A 
A 



sin[(^:5 - A)t^] ^ sm[{£s + A)t^] 



{£s - A) 



{Ss + A) 



sin^Ss-AYf] sm'[{Ss + Ay- 



i^s - A) 



{£s + A) 



(11.45) 
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I 1 2 

FIG. 2: The product \Fi{k,Q,p;t^)F2{k,Q,p;t^)\ vs. £3 in units of A for r^(A;) = 0.01 for 
^n{k)t^ <^ 1 for = 2,20 respectively. 



In the long time limit, in the first line we can replace 

7rS{Ss ± A) , (11.46) 



{Ss ± A) 

whereas the second term is subleading in the limit — )■ cxo; A — )■ 0, therefore I-F1-F2I — 
sin(A(p) t^)/(A(p) t^) |-Fi,2p, namely in eqn. flll.35p the interference term is of the 
same order as the direct terms for A(p) <C 1 but suppressed for A(p) ^ 1. 

We highlight that A{p) = n/tosc where tosc is the oscillation time scale, therefore the 
conclusion from the analysis above is that Fi ~ F2 when the smallest of TT^{k),t^ <^ tosc 
and the time uncertainty due to either the lifetime of the source or the disentanglement 
time scale leads to an energy uncertainty ~ VT,{k),2Ti /t^ that cannot distinguish between 
the mass eigenstates, and the resulting quantum state is similar to a "Pontecorvo" state 
namely a coherent superposition of mass eigenstates leading to oscillations with the familiar 
quantum mechanical survival probability. 

On the other hand if TT,{k),t^ > tosc the mass eigenstates separate in time as the energy 
uncertainty becomes smaller than the distance between the broadened "delta functions" , the 
coherence and interference between the mass eigenstates is suppressed by this separation and 
the disentangled state is not of the form of a "Pontecorvo state" . 

The energy uncertainty is determined by the smallest of the lifetime of the source and 
the disentanglement time scale. When this uncertainty is larger than the energy separation 
between neutrino mass eigenstates, the disentangled quantum states is the "Pontecorvo" 
state, but if the uncertainty is much smaller than the energy separation coherence and 
oscillations are suppressed. 

Therefore entanglement during a time scale of the order of or longer than the oscillation 
time scale suppress coherence and oscillations. 

These are some of the main results of this article. 

These results were obtained from the free evolution of the disentangled state, however, 
what is needed is an assessment of how these considerations are manifest in the energy 
spectrum of the charged leptons that are measured at the detector. 
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III. CHARGED LEPTON DETECTION: 



To obtain the transition probability to a final state with a charged lepton, which is 
the state finally detected, we must go to first order in the charged current interaction in 
flll.Sp . evolving the disentangled state \V^(t^)) up to time t^i at which the charged lepton is 
detected. In first order in the charged current interaction we find from flll.3ip 



iHn tj 



I •' 'A" 



to 



dt' / (fxUiiW{x,t')l{x,t')Pi{x,t') e'^"*" \V^{t^)). 



(III.l) 

The contributions that yield a W and a charged lepton / in the final state correspond to 
annihilating the neutrino states from |V^(t^)) and creating the final I. A straightforward 
calculation yields, 



|V,(t..))« 



W 9n)e~ 

sm{e) 
21^1 (p) 



cos(6) 

12 7^7r^^2[k,Q,P,tf,] G2[p,qi,tD,tf,] 



F,[k,Q,p,t,] Gi[p,qi,tD,t,] \ e-^^^*° \Wf:J 



(111.2) 



where 



p = k-Q ; k^=p-qr, Ed = Ew{k^) + Ei{qi) ; Es = E^{k) - E^{Q) (111.3) 

V" 2E^{k) 2E^{Q) 2Ew{k^) 2Ei{qi) ' 



n, 



and 



to 



Gi[p, qi, tn, g = / dt' e^(^^-^«(p))*' = ^KEn-o,ipmtn+t,)/2 



sin 


{ED-ni){tD- 






{Ed - a)/2 





(111.4) 

The functions F\^2 and G\^2 determine approximate energy conservation at the production 
(-^1,2) and detection (6*1,2) vertices respectively. 

The relevant observable is the energy distribution of the charged leptons measured in the 
detector, namely 

M{qi,tD) = «(V^(tD)|aJ(g1)az(g1)|V^(tB))« (III.5) 
using Qi{p) ~ E^{p) in the denominators, we find 

2 



N^,{q^,,tD) 



2 E,{p) 



cos\e) iFal^ |G2p+sin^(^) \Fi\^ \Gi\'^+- sm\2e) Re F2 F* G2 G* 

2 - - 

(111.6) 



Meiqe.to) 



2 E,{p) 



sm'^{29) 



I P |2 ir- |2 



+ ^\Fi\'\G,\' 



- Re 



F2 F* G2 Gl 



(III.7) 
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The term Re F2 G2 Gl describes the interference between the mass eigenstates that 

include the initial correlation in the entangled quantum state f lll.lSp . 

As a guide, we note that if \Fi \ = \F2\ and \Gi\ = \G2\ then one would find 



2 E,ip) 



gg-K V He 

2 E,{p) 



iFsHGar <! l-sin2(2^) sin^ 



I P |2 1/^ |2 

r 2I |tJ2| 



sin2(2^) sin^ 



(III.8) 



(III.9) 



where $ is the total phase in the product inside the real part in (]III.6|III.7|) . Clearly the 
terms inside the brackets in (]III.8IIIII.9I) are the disappearance V^i^^u^ and appearance V^^^ue 
probabilities associated with Pontecorvo states respectively. 

The general case is obtained by replacing F12 by (]II.16|) and Gi^2 by ( 1III.4p in the inter- 
ference term 



Re 



F2F*G2Gl 



(III.IO) 



Before we study this interference term in detail, let us consider again the direct terms 
|Fi^2p and |Gi,2p- The functions |-Fi,2|^ are sharply peaked at Es = ^i{p), therefore consider 
integrating either one of these functions with a density of states that is varies smoothly near 



\F,\^dEs ^ p{n,) 27cT^ik) 



(iii.ii) 



where we used eqn. (III.24P which imphes the identification flll.25p . Similarly for large tj;) — t^ 

/oo n 
p{En)m^dED^ p{n.:) / \G,\^dED-pm2'n{to~t^) (III.12) 
00 J 



hence just as in Fermi's Golden rule we identify 



\Gi\^ = 2'K {tD-t^)SiED-n,] 



(III.13) 



Furthermore for < 1 (eV)^ which is the putative mass range of sterile neutrinos to 
explain the short baseline anomalies, and typical neutrino energies E^, > 3 Mev with the 
lower range applying to reactor neutrinos, then A < 10"'' (eV) and A/E^ < 10"^'^. In 
typical neutrino experiments, the energy spectrum is measured with a finite resolution and 
"binned", namely integrated over an energy range determined by the resolution, however, 
such resolution is always much larger than 10~^ eV and in all experiments the relative error in 
energy resolution AE/E ^ 10~^^. The point is that the measurement resolution is nowhere 
near enough to discriminate an energy difference A{p) between the energy eigenstates in the 
"binning", and replacing 5{Es,d — Ey ± A) ^ 6{Es,d — Ey) is an excellent approximation. 
Therefore we can safely replace 



l^il 



2TxT^{k) 



5{Es-Ey{p)) =27rT^{k) 



5{Es) . 
(III.14) 
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= \G2? = 27r{tD - t^) S{Ed - E,{p)) = 27r{tD - t^) S{Sd) • 



(III.15) 



The approximations above rely on that Fi 2; (^1,2 ai'e distributions that are sharply peaked 
and they must be understood as being integrated with density of states, which experimentally 
are insensitive to the energy difference A(p) between the neutrinos. Under these approxi- 
mations, consider the product 



1-2 e*'^(p)*'' g-r^(fc)tM/2 cos (£^5^^) + e^*^^^')*'' e'^'^'^'^^^" 



(III.16) 



A{p) 



This function features two poles in the complex Es plane at Ss — ±(A(p) + ir^(A;)/2), 
being integrated with a density of states that is insensitive to A(p) in the narrow width 
approximation 

poo 

p{Ss)F2F* dSs p(0) / F,F* dSs 



f 

J —(. 



(III.17) 



the integral can be done and we find 



J — c 



F2F* dSs = 27T T^{k) 



1 + i 



2A(p) 

r,r(fc) 



1 + 



4A2(p) 



1 _ (,'^i^(p)ti. g-r^Wv 



(III.18) 



using these results for = we obtain 



p{Sn)G2GldSDC^p{0) (27r) e-^(^)(*-+*'') ^m)] 

A(p) 



(111.19) 



Therefore, under the assumption that the experimental "binning" is insensitive to the 
neutrino energy difference we can safely approximate the direct terms as. 



\Fi\' = 


\F2\' 


= 


\G2\' 



and for the interference terms 



(III.20) 



1 + in 









F2F* = 27iT^{k) 



G,Gl = 2.e--(^)(*-^V) sin[A(p)(t^-g] 

A(p) 



where we have introduced the ratio 

^_ 2A(p) _ 

T,{k) \2M^ro) E,{p) ' 

where Fq is the rest-frame decay width of the pion. 



(III.21) 



(III.22) 
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With these approximations, the number (density) of charged leptons ( muons ) mea- 
sured at disentanglement time (1II.18|) . is obtained by using the approximation (1II.25P and 

fli{p) ~ Ey{p), we find 



(III.23) 



9>VE^{k)E,,{Q)E,{p) 
and the number (density) of charged leptons measured at the detector is given by 

M^ito) = AfM dV^^w^ \ (cos^(^) + sin^(^)) {t " ^/.) + ^ sin^(2^)r[t,,, g ) (III.24) 



and 
where 



Kito) = Ar^(t,J dT.^we 2 sm\29) <j {t - t^) - Tito, 



(III.25) 



1 



1 



sin [A{p)(tD - 



cos[A(p)(tB + t^)] + n sm[A{p){tD + t^)] 



cos[A(p)(tB - t^)] + 7^ sin[A(p)(t^ - t^)] 



(III.26) 



and 



dr 



27T gHjEp - E,{p)) 
8VEw{k^)Ei{qi)E,{p) 



(III.27) 



is the differential charged lepton production rate from the reaction v ^ W I for a neutrino 
of energy E^ at the detector. 

These expressions become more familiar if we calculate the detection rate as is the usual 
procedure in S-matrix theory, we find the simpler results 



dM^jtp) 
dtp 

dXeito) 



dt 



D 



(111.28) 
(III.29) 



where the survival (disappearance) and appearance probabilities are 



cos[2A(p)tfl] + 7^sin[2 A(p)tD] - 



e-r.(fc)t^ cos[2 A(p) (t^ - t^)] + 7^ sin[2 A(p) {Id - t^)] 



(III.30) 
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■sm\29)\ 1 



1 



1 



1 + 7^2 1 



cos[2 A(p)tB] +7^sin[2 A(j))tz)] - 



cos [2 A(p) (tz) - t^)] + 7^ sin [2 A{p) {Id - t^)] 



(III.31) 



These expressions are in agreement with the previous discussion: when T^^ ^ 2 A(p) and 
sc ^ ^/i) the mass eigenstates cannot be discriminated during the hfetime of 
— )■ and the appearance and disappearance probabihties are given by the 



to > (or t 
the source, IZ 
usual result 



V^^^ito) = l-sin2(2^) sin' 



V^^eito) = sm\29) sin' 



D 



However, in the opposite limit 2A(p) ^ r,r(^), namely 7^ ^ 1 or to 



(III.32) 

the mass 



eigenstates are completely separated by the time evolution and the oscillation probabilities 
are suppressed. 

The origin of the discrepancy between between the probabilities (lIII.30|III.3ip and the 
familiar results given by (]III.32p is traced to the interference term Iie[F2F^G2Gl], the func- 
tions Fj,Gj are completely determined by the time evolution of the quantum state and 
describe the approximate energy conservation at the production and detection vertices. The 
functions Fi^2 describe the initial correlations in the entangled quantum state (111.13p . 

Therefore the physical reason behind the difference between the probabilities 
(lIII.30IIIII.3ip and (1III.32P is that the time scales associated with the lifetime of the source 
and the entanglement of the charged lepton define energy uncertainties which determine 
whether the mass eigenstates are separated during these time scales or not. Short lifetimes 
and disentanglement time scales (T^(fc) ; -C tosc) introduce large energy uncertainties and 
the mass eigenstates are "blurred" into a Pontecorvo state which yields the usual quantum 
mechanical result flIII.32p . In the opposite limit, long lifetime and disentanglement scales 



(T7r(fc) ; tfj_ ^ tosc) lead to small energy uncertainties and the mass eigenstates are separated 
leading to decoherence which is manifest in the expressions flIII.30IIIII.3ip in terms of TZ and 

It is important to highlight that the factorization in flIII.28|III.29p is a direct consequence 
of the fact that the binning or energy resolution in all current experiments cannot distinguish 
the energy difference between the mass eigenstates A(p) < 10"'^ eV and energy resolutions 
A{p) / Ey{p) < 10~^^ therefore the approximations (jlll. 1411111. ISP are amply justified. 

The interplay between the lifetime of the source and the disentanglement time scale and 
the suppression of the oscillatory component of the transition probabilities is more clearly 
exhibited in two simple cases: 

Case I: r7r(k)t^ ^ 1: This corresponds to a disentanglement time scale much larger 
than the lifetime of the source, in which case the energy uncertainty is determined by T.^. 
In this case 



sin'(2e)<^ 1 



l + 7^2 



cos [2 A(p) Id] + 7^ sin [2 A(p) to 



cos[2 A(p) to] + 7^ sin[2 A(p) to] 




(III.33) 



(III.34) 
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This expressions make clear that when A{p) ^ Tj^lk), namely for 7^ ^ 1 oscillations are 
suppressed in agreement with the analysis presented in the previous section. The suppression 
is a consequence of the separation of the mass eigenstates and the ensuing loss of coherence. 



Note that whereas the usual expressions (1111.320 valid for 7?. = are such that 



1 ; V^^e — )■ as ti) — )■ 0, for 7?. 7^ this is not the case, because the expressions (llll.33flll.34p 
only hold for t/j > ^ T^. 

Case II: r^(k)t^ < 1: 

In this case the disentanglement time scale is much shorter than the lifetime of the source, 
namely the source is nearly stationary during the time scale of disentanglement and we can 
simply approximate this case by taking TTj{k) — )■ 0. This approximation correctly describes 
the fact that the main energy uncertainty is determined by 27r/t^ as explained following 
equations (lll.42fll.43p . In this case we find 



jr. = e-*^s «m/2 gjf^j t^,/2 21 



sin 




- 




Es 







and invoking the same approximations as in the previous case we find 

2nt^6{Es-E,{p)) 



\Fi\' = \F2\' 



The interference term is given by 



2nglt^6{Es-E,{p)) 



Re 



F2F*G2Gl = cos(A(p)tB 



8VE^{k)E^{Q)E,{p) 

sin(A(p)t^) sin(A(p) {to-t^)) 



(III.35) 

(111.36) 
(111.37) 



A{p) A{p) 
X (27r) 6{Es - E,{p)) (27r) 6{Eo - E,{p)) , 



leading to the final results 



dt 



D 



where the survival (disappearance) and appearance probabilities are 



1 - -sin2(2^) 



■ sin (A(p)g 

1 i-7-^ COS 

A(p) 



t, 



■sm\29) 



sin(A(p)g 
1 — - cos 



A{p)t 



2Aip){tn-^) 



t. 



2A(p)(t^-|)j 



(111.38) 

(III.39) 
(111.40) 

(111.41) 
(111.42) 



As A{p)t^ — )• 0, namely when t^/tosc 0, the expressions for the probabilities become the 
familiar ones, but for A{p)t^ S> 1 the expressions above display two sources of suppression 
through entanglement: the prefactor sin (A(p) t^)/A(p) t^, and also a shortening of the 

effective baseline from L = cto to Lgff = c{tD — ^) ■ 
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IV. IMPLICATIONS FOR ACCELERATOR AND REACTOR EXPERIMENTS: 



The discussion of the previous sections hinges on two generation mixing, however, if 
sterile neutrinos are the correct explanation of the short-baseline anomalies then 



sm\29) -> 4|[/e4nf/; 



m4| 



(IV.l) 



It is convenient to write the probabilities in terms of the baseline L = ctj:, and introduce 
the disentanglement length of the muon = ct^, writing as usual 



A(p) to = 1.27 



L/m 



5m? 



eV' ^.(p)/MeV 



A(p)t^ = 1.27 



5m? 



eV' E.(p)/MeV 



and 



„5m2 MeV 
A(p) = 1.27— 2 



vr 



eV^ E^{p) Lo,c/m 
The disappearance and appearance probabilities are given by 



(IV.2) 



(1V.3) 



Va^JL) = l--sin2 20 <1 , 

^^^^ ^ 2 ^ M 1 + 7^2 1 _ g-Ld/L,(fc) 



cos [2 A{p) L]+n sin [2 A(p) L] 



COS [2 A(p) (L - Ld)] + n sin [2 A(p) (L - L,)] 



(IV.4) 



V.^eiL) = -sin2(2^)Jl L^- 1 

t^^^y I 2 \ l + 7^2l-e-^ 



COS [2 A (p) L] + 7^ sin [2 A (p) L] 



-La/LTy(k) 



cos[2 A(p) (L - Lrf)] + 7^ sin[2 A(p) (L - L,)] 



(1V.5) 



where L^^^k) = c/VT,{k) is the decay length^ of the parent particle (pion). 
For neutrinos produced from pion decay the ratio 71 (1111.22^ becomes 



n = 2n 



Lose {^Ey 



0.14 



E^{k) 
(eV)2 \E,{p) 



(IV.6) 



therefore, TZ 1 for 5m'^ ^ (eV)^ which is the range of masses for sterile neutrinos that 
could solve the short-baseline anomalies, and in the case of MiniBooNE 1 < E-j^jEy < 6. 

Remarkably, eqn. f lIV.4p is exactly the same as eqn. (23) in ref . [46|] where 7^ is equivalent 
to the quantity ^ and replaces the "pipeline" Id in this reference. Hence, the result of 
ref. [3] can be interpreted as disentangling the muon at the distance Id which is identified 
with the "pipeline" . 

Both LSND and MiniBooNE are designed with (L/m)/(E^/(MeV) ~ 1. 

At MiniBooNE a neutrino beam is obtained from pions that decay in a decay "pipe" 
~ 50 m long, and neutrinos go through ~ 500 m of "dirt" before reaching the detector 



^ For ETr{k) > 1 GeV the Lorentz factor 7 > 7 and we approximate /3 ~ 1. 
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determining a baseline L ~ 550 m with a peak energy in the neutrino spectrum at about 
600 — 1000 MeV. At MiniBooNE muons with ~ GeV energy are stopped at a distance ~ 4 m 
in the "dirt" thus La ~ 54 m^. 

Fig. ([2]) displays these probabilities for the set of parameters consistent with (one) sterile 
neutrino with Sm^ ~ leV^ ; sin^(20) = 0.2 and MiniBooNE baseline and range of neutrino 
energies. 




500 1000 1500 2000 2500 



FIG. 3: Appearance V^^e probability vs. Ei, (MeV) for MiniBooNE parameters. The solid line(s) 
correspond to pV.SP for E-,^(k) = 1,2 GeV (indistinguishable on the scale of the figure). The dashed 
line corresponds to the quantum mechanical probability V^^^e pil.32p . 



The figure shows that the appearance probability is suppressed as compared to the quan- 
tum mechanical result, the suppression being more pronounced at smaller energy where TZ is 
larger (see below). Although for 7^ 7^ the probabilities cannot be fit by the usual quantum 
mechanical result in the full energy range, a fit of the form 

V^^e{L) = sin2(2^e//) sin^ 

in a restricted energy range would lead to 

sin2(2^e//) < sm^{2e) ; Smlfj < Sm"^ (IV.8) 

as can be seen from the position of the maxima of the appearance probability: lower in 
amplitude (smaller mixing angle) and moving towards smaller energy (smaller dm?). For 
the case of MiniBooNE the fit is shown in figH] 

Because in this situation ~ L^^ik) <^ L decoherence from the source lifetime or 
entanglement does not lead to experimentally substantial corrections. 

Although not relevant for the MiniBooNE experiment, but as an illustrative example to 
display the effects of decoherence on the transition probabilities as a consequence of long 



^ The author is indebted to WiUiam C. Louis III for extensive correspondence clarifying these experimental 
aspects of the MiniBooNE experiment. 



4E. 



(IV.7) 
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Fit to the QM formula with eflFective parameters 



Ph 







L= 550 m 


/ \ 


E.-1 OeV 


La-54m 






sin^(2e) = 0.2 






Sm^=leV^ 






EflFective parameters for fit to QM 






sm^(2eeff) = 0.198 




1 


6m^etr=0.975 eV^ 

1,1 


500 


1000 1500 


2000 2500 






E(MeV) 



FIG. 4: Fit of P/,_^e vs. E,y (MeV) for MiniBooNE parameters and = IGeV. For sin2(26l) 
0.2, 5m^ = leV^ the fit yields sin2(26le//) = 0.198, 6mljj = 0.975 eV^ 



distance entanglement, we consider the case Ld ^ L^^ik), in which case the probabilities 
f HVaiV^I) simplify to 



V^^,{L) = 1 - ^sin2(2^)|l - ^-i^ ^cos[2 A(p)L] +7^sin[2A(p)L] 



1 + 



cos [2 A{p) L]+n sin [2 A{p) L] 



(IV.9) 
(IV.IO) 



These are displayed in fig. for L = 600m, ^m^ = leV^ sin2(2^) = 0.2, E^{k) 
1, 2 GeV. 





500 1000 1500 2000 2500 3000 



500 1000 1500 2000 2500 3000 



Ev(MeV) 



Ev(MeV) 



FIG. 5: Disappearance V^^fi and appearance V^^e probabilities vs. E^, (MeV) for L = 
600m,5m2 = leV^ , sin (20) = 0.2. The value E.,^ = refers to 7?. = 0, the usual quantum 
mechanical result for the the probabilities. 



Figure displays the appearance probability given by flIII.3ip . for the parameters E^, = 
2000 MeV ; ^m^ = leV^ ; L = 600 m ; sin2(2^) = 0.2 (sohd line) and the best fit to 
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the quantum mechanical probabihty fllV.Tp resulting in Smlff = 0.71 eV^ ; sm\29, 



0.185. Several aspects are clarified by this example: i) the suppression by lifetime and 
disentanglement effects leads to an underestimate of both 6m'^, sin^(2^), ii) the fit is reliable 
only within an intermediate energy range, much less reliable in the low energy region, iii) 
the ratio TZ implies that there are more parameters than the amplitude sin^(26') and the 
ratio L{m)/E^{MeV). 




Ev (MeV) 



FIG. 6: Disappearance Vfj,^^ for MiniBooNE parameters: the dashed line is the result with = 
2000 MeV ; Sm^ = leV^ ; sin (29) = 0.2, the solid line is a fit with the quantum mechanical 
probability (IRTtI) with 6mljf = 0.71 eV^ ; sin2(2(9e//) = 0.185. L = 600 m. 



Therefore, since the experimental data is always fit with the usual quantum mechanical 
formula, the values of sin^(26') ; from the fit actually correspond to sin^(26'e//) ; ^^l/f 
the above analysis leads to conclude that decoherence from the decay of the parent particle 
and the disentanglement of the charged lepton imply a larger value of the mixing angle and 
5m? from those extracted from the fit to the usual quantum mechanical probability. 

As shown above for the parameters of MiniBooNE, decoherence through lifetime and 
entanglement effects yield very small corrections, however the principal and fundamental 
observation remains, namely lifetime or disentanglement time scales similar to or larger than 
the oscillation time scale lead to decoherence and suppression of the appearance probabilities. 
A quantum mechanical fit yield effective values 6'e// , ^^Ijf which are smaller than the actual 
values. 

In our analysis we have assumed that the entangled quantum state arises from the two 
body decay of a parent particle, (here considered to be the pion), however at LSND the 
(anti) neutrino beam is produced by the three body decay of a muon at rest, whereas at 
reactors the (anti) neutrinos are produced via the /3 decay of long-lived unstable nuclei 



235^ ^ 238^ ^ 239p^ ^ 2^^Pu[37|, [SSj. Although the actual calculation presented in the previous 
section for the exact entangled state does not directly apply to the description of the quantum 
states of neutrinos produced at LSND and of reactor experiments, in absence of a more 
detailed understanding of the entangled quantum state resulting from the three body nuclear 
P decay, we will use the result (1II.13P with the caveat of possible corrections arising from 
three body phase space effects. 

At LSND muon antineutrinos are produced from vr"'" — )■ /i^ i/^ followed by e^VeV^ 
where most of the muons decay at rest. The resulting beam attains the maximum 
energy at the Michel end point 52.8 MeV and the liquid scintillator detector is located 
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about L = 30 m from the neutrino source. Since L <^ cr^ ~ 660 m for LSND it follows 



that r^t^ <C 1 (the relevant decay width now is the muon's as the parent particle) and 
this corresponds to case II (nearly stationary case) with Ld the disentanglement length. 
The same limit applies to reactor experiments where neutrinos are produced from nuclear 
P decay of long-lived radiaoactive nuclei, therefore for LSND and reactor experiments the 
disappearance and appearance probabilities are given by, 



1 - -sm^{2e) 



■ sin (A(p) Lrf) 

1 W N r COS 

A(p) Ld 



2A(p)L 1 



2L 



(IV.ll) 



sin (A(p) Ld) 

1 A / X -r COS 

A(p) Ld 



2A(p)L 1 



Ld_ 
2L 



(IV.12) 



In ref. [2j] it was also recognized that the muon lifetime does not affect the transition proba- 
bilities at LSND, however, the effect of disentanglement has not been previously recognized. 
In LSND, the detector is at L = 30 m from the neutrino source and is shielded by the 



equivalent of 9 m of steelj48j which then should be taken as a figure of merit for Ld ^ 20 m. 
At reactor experiments a figure of merit could be the size of the reactor core, at CHOOZj49| 
it is approximately ~ 4m with a baseline 100m ^ L < 1km, although, quite likely these 
figures of merit for Ld overestimate the disentanglement length scale both in LSND and in 
reactor experiments. Unlike the case of MinibooNE where the suppression factor is mainly 
determined by the pion decay length, at LSND and reactor experiments the disentanglement 
scale Ld is less certain. 

Thus we take Ld as a parameter and study the disappearance and appearance probabilities 
within the range Om < Ld < 15m to illustrate the consequences of decoherence from 
entanglement and to extract the main conclusions. These are displayed in fig. ([7]). 




FIG. 7: Disappearance V^-^fi and appearance P^-^e probabilities vs. Ei, (MeV) for LSND and 
reactor parameters. The value Ld = refers to the usual quantum mechanical result for the the 
probabilities. 



These figures reveal a situation very similar to that analyzed above for MiniBooNE. 
Larger disentanglement lengths Ld lead to a larger suppression of the appearance probabil- 
ity. Similarly, a fit to the experimental data with the usual quantum mechanical appearance 
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probability results in an underestimate of both sin^(2^) and Sm'^ for the same reasons ana- 
lyzed above. 



V. COMMENTS ON WAVE PACKETS: 



The study in this article was restricted to plane waves to exhibit the main results and 



conclusions in the clearest possible setting. As has been argued in the literature [16|, 
20|, H HI, M wave packet localization is an important ingredient in the description of 



neutrino oscillations. The localization length both of the production and detection regions 
define momentum uncertainties that are important in the conceptual understanding of the 
interference phenomena. 

Furthermore, in our calculation the disentanglement t^ and detection tjj times are sharp, 
this is a consequence of calculating the transition matrix elements in finite time intervals, 
however, the wave packet treatment smears these times over the time scale during which the 
wave packet overlaps with the detectors which is the appropriate physical description of the 
detection events. 

The analysis in[l^, 17, 20, 23, 25, 26, [s^ (typically with Gaussian wave packets) clar- 
ifies that neutrino wave-packets evolve semiclassically, the center moves as the front of a 
plane wave with the group velocity and is modulated by a Gaussian envelope which spreads 
through dispersion. Wave packets associated with the different mass eigenstates separate as 
they evolve with slightly different group velocities and when their separation becomes of the 
order of or larger than the width of the wave packet the overlap vanishes and oscillations are 
suppressed, typically exponentially in the ratio L^/L^^^^ where Lr„h — a El/5rn? and a is the 
spatial localization scale of the wave packet. As discussed in |5l| the wave packet descrip- 
tion also features another source of decoherence in the localization term, which suppresses 
coherence when a > Lose- 

However, it should be clear from the discussion and results presented above, that energy 
uncertainties from the width of the parent particle, disentanglement time scales, finite time 
intervals between production and detection and experimental measurements are sufficient 
to guarantee interference and oscillations. Entanglement over long distances and time scales 
introduces decoherence in a quantifiable manner. Introducing wave packets will modify 
the results only quantitatively but by no means fundamentally: a wave packet is a linear 
superposition of plane waves and the analysis for each plane wave described above can be 
generalized to such superposition. One aspect that relies on a wave-packet description is the 
detection: the total number of events is obtained by the event rates multiplied by the total 
time that the wave packet takes to pass through the detector. For ultrarelativistic neutrinos 
this is of order a/c since spreading through dispersion can be neglected on short baselines, 
therefore the total number of events is given by the rates (llll.39pil.40p multiplied by a/c, 
obviously this will not change the distortion of the spectrum determined by the oscillations 
in the appearance and disappearance probabilities. Another correction is the geometric 



fiux factor which again for short baselines can be neglected. As found in refs. [18|, |2J, |34 



including the lifetime of the source in the wave-packet evolution introduces another length 
scale (the decay length of the parent particle) which competes with the localization length 
of the wave-packet. As discussed above, wave packet localization will not affect oscillations 
unless the wave-packets corresponding to the different mass eigenstates begin to separate. 
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For Sm^ ^ 1 eV^ ; Ei, ~ MeV and L ~ 100 m the criterion for separation over the basehne 
would require a locahzation length a < 1 A, this estimate is much larger than the nuclear 
radius for unstable nuclei, thus decoherence via the separation of the wave-packets of mass 
eigenstates may be another source of decoherence if the localization length scale of the wave 
packets is of nuclear dimensions. 

Thus we conclude that the results obtained with the plane wave analysis will apply vis 
a vis to the case of wave packets, unless the source of decoherence associated with the 
separation of wave packets of mass eigenstates introduces enough decoherence as to dwarf 
the effects discussed here. On the short baseline experiments considered here this would 
require localization lengths < 10^^° m for reactor experiments and < 10^^^ m for accelerator 
experiments. 

Strengthening these arguments require (and warrant) a full study of the complete de- 
scription of disentanglement and lifetime effects in a wave packet formulation. Of particular 
importance is whether for Sm? ~ eV^ ; ~ few MeV wave packet localization on nuclear 
scales can be a source of decoherence in reactor experiments. The results of this study will 



be presented elsewhere [47 



Wave packets vs. disentanglement: 

Decoherence through lifetime and disentanglement is fundamentally and conceptually 
different from decoherence in the wave packet formulation. Neutrino wavepackets manifestly 
describe single particle states that are spatially localized, the spatial localization introduces 
uncertainty in the momentum, and in this formulation decoherence is a consequence of the 
separation in space of the wave packets associated with the different mass eigenstates. As 
explained in ref. j5l| there are two sources of decoherence: one resulting from the separation 
of the wave packets of different mass eigenstates through their differe nt g roup velocity, and 
another determined by a localization term (see eqn. (8.114) in ref. [l2|) which results in 
decoherence for a ^ Lose- 

Entanglement, on the other hand, refers to the fact that the quantum state that results 
from the decay of the parent particle is a correlated many particle state, the correlation be- 
tween the charged lepton and the neutrino(s) is manifest in the coefficients Ci{k, q,p, t) in the 
quantum state f lII.13|) . These coefficients are time dependent and describe the approximate 
conservation of energy at the production vertex. A single particle neutrino state is obtained 
by projection of the charged lepton state, this projection is the quantum mechanical mani- 
festation of the observation, absorption or decay of the charge lepton and disentangles the 
(two body) quantum state at a time scale t^. These correlations are precisely the origin of 
the terms Fi 2 which enter in the interference term (llll.lOp and are, therefore, the origin of 
the difference with the familiar quantum mechanical result. In this description the lifetime 
of the source and determine energy uncertainties as explained in the previous sections. 
Decoherence ensues when the energy uncertainty is much smaller than the energy separation 
between the mass eigenstates. This source of decoherence is obviously independent of the 
spatial localization of the quantum state and is present even for plane waves, unlike wave 
packet decoherence. Although decoherence in the wave packet and disentanglement formu- 
lations are physically and conceptually different, they are indeed complementary and both 
will be present in a complete wave packet description of neutrino oscillations. For example 
as discussed in ref. [5l[ if a neutrino wavepacket produced by the decay of a parent particle 
of width r is assigned a localization length 1 /F then the condition for decoherence from the 
localization term a ~ Lose becomes equivalent to F ~ A(p) which is the condition which 
results from the disentanglement analysis in the case when the lifetime is shorter than the 
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disentanglement time scale. However, obviously this cannot be the case for reactor neutrinos 
since the lifetime of the parent particle is thousands of years and the relevant scale is the 
disentanglement length scale as discussed above. 

VI. CONCLUSIONS AND FURTHER QUESTIONS: 

Accumulating evidence for anomalies in short-baseline experiments pointing towards a 
change in the current paradigm of neutrino oscillations resulting from the mixing among 
three active species, will likely motivate further accelerator and reactor short baseline ex- 
periments. The firm assessment of new "sterile" neutrinos as possible explanations of the 
data warrant a deeper understanding of quantum coherence that determine the appearance 
and disappearance probabilities. 

The realization that the neutrino states produced in charged current interaction vertices 
are quantum entangled states of the neutrino and its flavor charged lepton partner call for 
a re-examination of the usual quantum mechanical description of neutrino oscillations as 
simple two level systems (for two neutrinos mixing). The measurement, absorption or decay 
of the charged lepton leads to the disentanglement of the quantum state, but the resulting 
neutrino state features the correlations from the prior entanglement. 

The disentanglement of this correlated quantum state is a necessary condition for coher- 
ence between the mass eigenstates leading to oscillations, entanglement over long time scales 
project out energy eigenstates preventing oscillations. The usual "Pontecorvo" (quantum 
mechanical states) emerge if the disentanglement time scale is much smaller than the os- 
cillation scale. This is a consequence of the time-energy uncertainty: for disentanglement 
time scales shorter than the oscillation time, the uncertainty in energy cannot discriminate 
between the different mass eigenstates, the longer the entanglement time scale the smaller 
the energy uncertainty and the mass eigenstates become sharply defined in the correlated 
state leading to a suppression of the oscillation probability. 

In this article we find that both the entanglement with the charged lepton and the lifetime 
of the source that produces the neutrino beam lead to a suppression of the appearance 
probabilities. The relevant dimensionless parameter that quantifies decoherence by both 
effects is the ratio tiLs/ Lose where Lg is the smaller between the decay length of the parent 
particle (source) and the disentanglement length scale. 

We obtain the corrections to the disappearance and appearance probabilities both from 
entanglement and lifetime effects in a model which captures in a clear and reliably manner 
the main features of the production, evolution and detection of mixed states. 

For MiniBooNE, the most important source of suppression is the decay length of the pions 
that produce the neutrino beam which is of the same order as the disentanglement length 
for the muons, whereas at LSND and reactor experiments, the disentanglement distance is 
the relevant scale that determines the suppression, for LSND this is because neutrinos are 
produced by muons decaying at rest while in reactor experiments neutrinos are produced 
via (3 decay of long lived radioactive sources, in both cases the disentanglement time scale 
is shorter than the lifetime of the source. 

Short baseline experiments imply small Lose therefore the impact of disentanglement and 
source lifetime is larger in these experiments. The suppressions of the oscillation probabilities 
are more pronounced at lower energies and are more dramatic for 5m^ ~ 1 eV^ which is the 
mass range for "sterile" neutrinos proposed as possible explanations of the short-baseline 
anomalies. 
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Our main results are the general disappearance and appearance probabilities given by 
eqns. f lIV.4|IV.5p . These simplify to equations (1IV.11PV.12|) when the disentanglement time 
scale is much shorter than the lifetime of the source, this is the case in reactor experiments 
(neutrinos at reactors are produced by (3 decay of long lived radioactive nuclei) and at LSND. 
The determination of the scale is cleaner in accelerator experiments where the neutrino 
beam is produced by pion decay (either at rest or in flight). However, for MiniBooNE the 
corrections are relatively small because the disentanglement length scale is of the order of 
the pion decay length and both are much smaller than the baseline. In reactor experiments 
Ld is more difficult to establish, a figure of merit is the size of the reactor core, but this 
estimate is probably too simplistic and overestimates the disentanglement length. 

While the experimental impact of the corrections in current experiments is relatively 
small, this work suggests that in the analysis of the data, the issue of disentanglement length 
scale must be addressed for a consistent interpretation of the results. An important corollary 
of our results is that fitting the experimental data with the usual quantum mechanical 
expressions for appearance and disappearance probabilities underestimates both sin^(26') 
and 5m^, furthermore this fit to the data differs substantially at low neutrino energy from 
the correct expression for the probabilities that include both the lifetime and disentanglement 
suppression, since the suppression is larger at smaller energies (shorter Lose)- 

An aspect that remains to be explored further is the description of neutrino propagation 
in terms of wave packets: the source and detector are spatially localized, in particular the 
localization of the source entails that the neutrinos are produced in entangled wave packets, 
the disentanglement of the charged lepton brings in another localization scale (at which 
the charged lepton is measured, absorbed or decays) which also infiuences the disentangled 
neutrino state. Wave packet localization also introduces yet another decoherence length 
scale Lcoh oc aE'^/Sm^ where a is the spatial localization scale of the wave packet. For 
Sm"^ ~ 1 eV^ sterile neutrinos in reactor experiments it is possible that Lcoh ^ Lose which 
would result in yet another source of decoherence and suppression of oscillations. These 
aspects are currently being studied and will be reported in a forthcoming studyj47|. 

Finally, it is worth commenting that quantum entanglement is also ubiquitous in B-meson 
oscillations, where the process of "fiavor" tagging actually disentangles the entangled B^ — B^ 

state produced by T(4s) decayjH, 52, Hsf, and quantum entanglement of the C = —\ Bg Bg 
pair produced in the decay of the T(5s) has been invoked for a determination of the width 
difference [s^] . Thus neutrino mixing is yet another fascinating manifestation of quantum 
entanglement in a system that maintains macroscopic quantum coherence over scales of 
kilometers. Fascinating examples of quantum entanglement on macroscopic scales are also 
emerging in other unlikely systems: photosynthesis in light harvesting complexes 551] and per- 
haps most surprising and provocative, as a possible explanation of the "avian compass" |56| . 
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Appendix A: The Wigner-Weisskopf Method 



For completeness we give a detailed presentation of the field theoretical version of the 
Wigner-Weisskopf approximation as it is not widely available in the literature. 

Consider a system whose Hamiltonian H = Hq+Hj where Hq is the free field Hamiltonian 
and Hj the interaction. The time evolution of states in the interaction picture of Hq is given 
by 

4t\^it))i = Hi{t)\^{t))j, (A.l) 
where the interaction Hamiltonian in the interaction picture is 

Hi{t) = e'^^'Hie^'^'^' (A.2) 

This has the formal solution 

\^{t))j = U{t,to)\^{to))i (A.3) 
where the time evolution operator in the interaction picture f/(t, to) obeys 

i^U{t,to) = Hi{t)U{t,t,). (A.4) 
at 



Now we can expand 



\'^{t))i = Y.^n{t)\n) (A.5) 



where \n) form a complete set of orthonormal eigenstates of in the quantum field theory 
case these are many-particle Fock states. From eg. (lA.ip one finds the exact equation of 
motion for the coefficients Cn{t), namely 

Cn{t) = -tY,Cm{t){n\Hj{t)\m) . (A.6) 

m 

Although this equation is exact, it generates an infinite hierarchy of simultaneous equa- 
tions when the Hilbert space of states spanned by {|^)} is infinite dimensional. However, 
this hierarchy can be truncated by considering the transition between states connected by 
the interaction Hamiltonian at a given order in Hj. Thus consider the situation depicted in 
figure IE] where one state, \A), couples to a set of states which couple back to \A) via 

Hi. 




FIG. 8: Transitions 1^4) o |k) in first order in Hj. 
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Under these circumstances, we have 



CA{t) = -iY,{A\Hi{t)\K)C^{t) 

K 

C.{t) = -tCAit){f,\Hj{t)\A) 



(A.7) 
(A.8) 



where the sum over k is over all the intermediate states coupled to \A) via Hj. 

Consider the initial value problem in which at time t = the state of the system = 
0)) = Ca{0)\A), namely 

C,{t = 0) = 0. (A.9) 
We can solve eg. flA.Sp and then use the solution in eq. flA.7p to find 



CAit) 



{K\Hj{t')\A)CA{t')dt' 
j:A{t-t')CA{t') dt' 



where, using ( ]A.2|) we find 



(A.IO) 
(A.ll) 

(A.12) 



This integro-differential equation with memory yields a non-perturbative solution for the 
time evolution of the amplitudes and probabilities. Inserting the solution for Cyi(t) into 
eq. flA.lOp one obtains the time evolution of amplitudes C,^{t) from which we can compute 
the time dependent probability to populate the state \k), \Cf^(t)\'^. This is the essence of the 
Weisskopf-WignerjSTII non-perturbative method ubiquitous in quantum optics 58| and the 
decay formalism of Kq — Kq mixing (59l|. 

The hermiticity of the interaction Hamiltonian Hj, together with the initial conditions 
in eqs. (]A.9P yields the unitarity condition 



1 . 



(A.13) 



Equation (lA.lip can be solved exactly via Laplace transform (60j, however, in weak cou- 
pling, the time evolution of CA{t) determined by eq.( ]A.ll| is slow in the sense that the time 
scale is determined by a weak coupling kernel S oc Hj. This allows us to use a Markovian 
approximation in terms of a consistent expansion in derivatives of Ca 60||. Define 



Wo{t,t')= f EA{t-t")dt" 

Jo 



so that 



d 



^A{t-t') = —Wo{t,t'), Wo{t,0) = 0. 
Integrating by parts in eq.l lA.ll]) we obtain 



(A.14) 



(A.15) 



EA{t - t') CaH') dt' = W^it, t) CaH) - Wo{t, t') ^CA{t') dt'. (A.16) 
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The second term on the right hand side is formally of fourth order in Hj and we see how a 
systematic approximation scheme can be developed. Setting 



Wi{t,t')= [ Wo{t,t")dt", Wi{t,0) = 
Jo 



and integrating by parts again, we find 



d 



Wo{t,t') -^CA{t') dt' = W^{t,t) CA{t) + 



leading to 



t') CA{t') dt' = W,{t, t) CA{t) - H^i(t, t) CA{t) + 



(A.17) 



(A.18) 



(A.19) 



This process can be implemented systematically resulting in higher order differential 
equations. Up to leading order in this Markovian approximation the equation eq. flA.lip 
becomes 

CA{t) [1 - Wi{t,t)] + W^{t,t)CA{t) = (A.20) 

with the result 

-iWo{t,t) 



l-W,{t,t) 



-iWoit.t) [l + Wi{t,t) + 



(A.21) 



To leading order in the interaction {0{Hj)) we keep S{t) = —i W^lt, t). Note that in general 
S{t) is complex. In the long time limit and using the representation flA.12p we find 



where 



(Ea - E,) 

is the energy shift in agreement with second order perturbation theory, and 



(A.22) 



(A.23) 



(A.24) 



this result for the width is in agreement with Fermi's Golden rule. Finally, in the Markovian 
approximation the Wigner-Weisskopf method yields 



CA{t)=CA{^)e 



(A.25) 



This solution agrees with the exact solution via Laplace transform 60|^. Inserting this result 
into equation flA.lOl) we find 



Jo 



-Ca{0) {^\Hj{0)\A) 



1 _ ^-i{E';^-E^~i^)t 



(A.26) 



Here we neglect wave function renormalization as it is not relevant for the discussion. 
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where = Ea + ^Ea is the renormahzed energy. The Schroedinger picture state |^(i))5 = 
e~*^°*|^(i))/ is finally given by 



E\-E^- 



iE^t 



For t ta = 1 /Ta the asymptotic state becomes 



\nt:>TA))s--CA{Q)Y. 



{K\HimA)er 



iE^t 



Fa 



k) 
(A.27) 

(A.28) 
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